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We show that the limiting distributions of the coefficients of the q- 
Catalan numbers and the generalized g-Catalan numbers are normal. 
Despite the fact that these coefficients are not unimodal for small n, 
we conjecture that for sufficiently large n, the coefficients are unimodal 
and even log-concave except for a few terms of the head and tail. 
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1 Introduction 

The main objective of this paper is to show that the limiting distribution of 
the coefficients of the g-Catalan numbers is normal. The Catalan numbers 



have many combinatorial interpretations, see Stanley [9]. The usual g-analog 
of the Catalan numbers is given by 
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C n (q) : 



1 ~2n 



(1.1) 



[n + 1] n 



where [n] = 1 + q + q 2 + ■ ■ ■ + q n 1 , and 

[k]\[n-k]V 

There are also other types of g-analogs of the Catalan numbers, see, for 
example, Andrews |2j, Gessel and Stanton [I], Krattenthaler [5]. 

We also consider the limiting distribution of the coefficients of the quotient 
of two products of g-numbers, which includes the result for the g-Catalan 
numbers as a special case. We conclude this paper with two conjectures 
on the unimodality and log-concavity for almost all the coefficients of the 
g-Catalan numbers and the generalized g-Catalan numbers provided that n 
is sufficiently large. 



2 The Limiting Distribution 

In this section, we use the moment generating function technique to obtain 
the limiting distribution of the coefficients of the g-Catalan numbers. We 
introduce the random variable £ n corresponding to the probability generating 
function 

0n(<?) = C n {q)/C n . 

As far as the computations are concerned, we will not need the following 
combinatorial interpretation of C n (q). However, for completeness and for 
the sake of presentation, we would mention that £ ra reflects the distribution 
of the major indices of Catalan words of length 2n, see, for example, [3]. We 
write 

Cn(q) = ^m n (k)q k . 
The following lemma is concerned with the expectation and variance of 

Lemma 2.1 We have 

, , n(n — 1) , , n(n — l)(n + 1) 
E(£n) = K 2 J and Var(£„) = ^ ^ >-. (2.1) 



n 
k 



2 



Proof. By the definition of C n (q), it is easy to check the following symmetry 
property of m n {k): 

n^n{k) = m n (n(n — 1) — k). 

Hence 

= 

Let 

n— 1 n—1 

F = F{q) = Y[{l + q + --- + q n+i ) and G = G(q) = JJ(1 + g + • • • + <f)- 

i=l i=l 

It is easily verified that C n (g) = F/G. Since 

/F^ _ FCT _ 2G"F' 2G' 2 F \ 
M?J l^i - ^ G G 2 + 



we obtain 



^n(n-l)(3n 2 -n-4)-C n , 



WU , ^ ^2 1 



Var(^) = " ' ' g - L + E(£ n ) - E^ n ) 2 = -n(n - l)(n + 1). 
C n 6 

This completes the proof. 



Lemma 2.2 VF/ien n — > oo, we /iawe 

fe=2 v ' i=2 

uniformly for t from any bounded set, where Bj 's are the Bernoulli numbers 
and a 2 is the variance of £ n as given in \2.1\) . 



Proof. The second summation can be expanded as follows: 

t «» + - **) -EE (?) = E (?) ( E 

j=2 i=2 j=l ^ J ' j=l ^ J ' \ i=2 



For k > 1, the second factor in the preceding summation is bounded by the 
following integral: 



i=2 

Consequently, 



« rn+l i , i\2fc-j+l _ 1 



2k-j + l 



((n + *) 2fe - * 2fe ) < 2 2fe (n + l) 2fc+1 < 8 2fe n 2fc+1 . 

i=2 

Since a 2 = > ^ when ra is sufficiently large, we have 

i=2 

for large n and k > 1. Thus 



2ffc,fc>3 i=2 



< n 



-1/3 



^ \ B 2k\ 



i2k 



2ffc,fc>3 



2A:(2fc) 



.g 4 2fc n -fc/3 



-1/3 



(64tn 



-iN2fc 



2A:(2fc)! 



2ffc,A:>3 

In view of the following asymptotic expansion of the Bernoulli numbers 

2(2n)\ 



B 



2n\ 



(27T) 2 "' 



the convergent radius R of the series 2J |-^2fc| 2 k(2k)\ e( l ua ^ s 2vr. Since £ is 

2\k,k 3 

from a bounded set, when n is large enough, the series 



l 52fc l 



-i\2k 



2\k,k>3 



(64tn 
2k(2k)\ 
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converges. Moreover, it is evident that QAtn a < 1, we can bound the above 
summation by the constant 



2Jfc(2Jb)!' 



Mi = V \B 2k \— - 

2ffc,fc>3 

Similarly, it can be deduced that 

2|fc,fc>2 v y i=2 

where M 2 = B 2 h 2 k(2k)i * s a cons t an t- Hence 

2|fc,fc>2 

fc=2 v ' i=2 

which tends to zero as n —>■ oo. This completes the proof. 



Theorem 2.3 When n — > oo ; £/ie random variable 

= — 

Var(£„) 5 

/ias i/ie standard normal distribution. 



Proof. Let M n (q) denote the moment generating function of £ n . Then we 
have M n (q) = <p n (e q ), see Sachkov [7]. Hence 

n + 1 l-e 9 A 1 - e(n+i)<? 

In/ i=l 

i JL 1 _ e (n+i)9 



TT— TT — 

i=2 i=2 

JL (l_ e (n+0ff)/( n + i ) 

11 (l-e*)/i 
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I " } n / e {n+i)q/2 e -(n+i)q/2\/n±i 

<^1> ^ 2 E (( n + 05 - «f) j II (e V2_ e -V2)/i 2 

eXP l 2 /II sinh(^/2)/| 
Recalling the following relation on the Bernoulli numbers [6] 

sinh(x/2)\ x 



we find that 



~Vf>A £(<»+••>" -•■")■ 



2 ^ 2fc(2fc) 



i=2 



Setting q = t/a, where a is the standard deviation of £ n as given in Theorem 
12.11 we are led to the expansion 

/ \ oo 2k n 

lnM„( t/ .) = ^ + £ S «" + <>" " ^ ' 

fc=l v 7 i=2 

Applying Lemma 12.21 we have, when n — > oo, 

00 4.2k n 

uniformly for t from any bounded set. Finally, 

n(n — l)t 



lim M n (t/a) exp 



2(7 



/.2A: n 



6 



t 2 

lim exp { B 2 —— £ ((n + if - i 2 ) 



•2(2)!^ i 2 



which coincides with the moment generating function of the standard normal 
distribution. Employing Curtiss' theorem [7], we reach the conclusion that 
r] n has the standard normal distribution when n approaches infinity. I 



3 A General Setting 

In this section, we will determine the limiting distribution of the coefficients 
of a quotient of products of g-numbers and will give two special cases. 

Theorem 3.1 Let ax, a%, a^, . . . and b\, b 2 , 63, . . . be two sequences of positive 
numbers, and let 

, ,n V- ,,\ k (1 ~ q ai )(l ~ q a2 ) ■ ■ ■ (I ~ q an ) 

<Pn(x) = > Pn(k)x k = ^-r^) -7— r 7 Vv- 

v v ; (l-g 6i )(i-g b2 )---(i-g 6 ' 1 ) 

Suppose that £ ra is the random variable corresponding to the generating func- 
tion <p n {x), that is, 

k 

Then £ n is normally distributed as n — > oo ; i/ and only if for k > 1 
^ t 2k f n 



8=1 



Proof. The expectation of £ n is easy to compute, as given below: 

1 n 

E(Cn) = M X )'<t=l = 9 E ^ ~ ^ • 



2 

1=1 



Proceeding analogously as in the proof of Theorem 12. 1} we find 

1 - 

o 2 = Var(£ n ) = - £ K 2 - b i) • 
i=i 

Hence, 



By the same procedure as in the proof of Theorem 12.31 we obtain 
lim M n {t/a)e W \ l -j^{af -bf) 



2 

i=i 



t 2k 



e«VM^ exp ^ B2t _J_,g (of _ a 



It follows that the limiting distribution of p n {k) is normal if and only if 



i2fc / n 



E B *2k(2k^ (E W - J-" as „ oo, (3.2) 

for £ from any bounded set. By virtue of the variance formula (13. ip . the 
condition (I3.2p is equivalent to 



\i=l 

for i from any bounded set. Thus f)3.2p is verified. This completes the proof. 



Corollary 3.2 Let p n (k) be given as in the above theorem. Suppose that for 
k > 2, there exist constants a > 0, j3 < and 7 < such that 



E W* - bf) 
=1 

EK 2 -& 2 )) 
i=i / 



^ £ < ^(an") J *, (3.4) 



fort from any bounded set. Then the limiting distribution ofp n (k) is normal. 



Proof. Note that the convergent radius R of the series 



E \ B2k \ 



x 



2k 



2ffc,fc>3 



2k(2k)\ 



is 2tt. If ( 13. 4p holds for k > 1, then for t from any bounded set, and for 
sufficiently large n, we have 



t 2kJ2( a f-bf) /a 2k 



i=i 



< n<{tan p f k 



where tan 13 < 2n. It is clear that n f — >• since 7 < 0. I 

If we choose a = 32\/3/3, 2/3 = 7 = — 1 ; Theorem 13.21 contains Theorem 
12.31 as a special case. We now give two more examples. One is the following 
g-analog of the Catalan numbers 

Cn(q) 



2n 



[2n] [n - 1 

which are symmetric and unimodal, see Stanley []. 



Using Theorem I3.1[ we reach the following assertion. 



Corollary 3.3 The distribution of the coefficients in c n (q) is asymptotically 
normal. 



Proof. First, we write c n (q) in the following form: 

n 

n(l-g n+i_1 ) 

i=3 

ra-1 ' 

(1 - g ) n (i - ?*) 

i=3 

Set a\ — d2 = 1, cbi — n + i — 1, 3 < i < n, and b\ = b 2 = 1, 63 = 1, hi = 
i — 1, 4 < i < n. Then we have 

n n 
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n-1 
i=3 



and 



n \ I n—1 



= (n + 2) 2 -l + ^((n + ,) 2 -. 2 ) 

,i=l / \ i=3 

= (n-l) fe (n + l) fc (2n-3) fc . 

By the same arguments as in the proof of Lemma l2.2l we may set a = 32^/3 
and 2/3 = 7 = — | such that the condition (13.41) is satisfied. Therefore, 
Theorem 13. II implies the limiting distribution of the coefficients of c n (q). I 

The m- Catalan numbers are defined by 

1 / mn 



n,m / -i \ i I 

(m — l)n + 1 \ n 



for n > 1. Accordingly, the generalized g-Catalan numbers are given by 



mn 
n 



[(m- l)n+ 1] 
Theorem 13. II has the following consequence. 



Corollary 3.4 The coefficients of the generalized q- Catalan numbers C n ^ m (q) 
are normally distributed when n — > oo. 



Proof. First, express C„ )m (g) as follows 



(m—l)n+i 



i=2 



Set ai = 1, Oj = (m — l)n + z, 2 < i < n, and &i = 1, 6j = i, 2 < % < n. 
Then we have 



n n n 2k ✓ v 

E(«f-if) = EW'-?) = EE ; 

i=l i=2 i=2 jf=l " ' 



((m — l)n) 2fc J ' 
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The same argument as in the proof of Lemma \2.2\ yields the following bound 
£(af-&f)<8 2fc ((m-l) 



s2k+l 



i=i 



Now 



EK 2 - & ') = £(((m-l)n + i) 2 -i 2 ) 

,i=l / \i=2 

> (m-l) 2k n 2k (n-l) k 

> (m - l) 2fe+1 n 3 V(2m) fc . 



It follows that 



V n (a 2k - b 2k ) 



(£ILi(af-& 2 )) 



2\\ k 



Again, by the same arguments as in the proof of Lemma I2.2[ we may set 
a = 8y / 2m and 2(3 = 7 = — | such that the condition (13. 4p holds. Finally, 
we may use Theorem 13.11 to get the desired distribution. I 



4 Open Problems 

While the g-Catalan numbers are not unimodal for small n, see Stanley [8], 
the limiting distribution suggests that the coefficients are almost unimodal 
in certain sense for sufficiently large n. Obviously, the first and the last term 
should not be taken into account otherwise one can never expect to have 
unimodality. In fact, an easy computation indicates that C n (q) are unimodal 
for n > 16. 



Conjecture 4.1 The sequence {m n (l), . . . ,m n (n(n — 1) — 1)} is unimodal 
when n is sufficiently large. 



When n > 70, numerical evidence suggestive of a stronger conjecture: 
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Conjecture 4.2 There exists an integer t such that when n is sufficiently 
large, the sequence {m n (t), . . . ,m n (n(n — 1) — t)} is log-concave, namely, 

(m n (k)f > m n (k + l)m n (k - 1) 

for t + 1 < k < n(n — 2) — t — 1. Moreover, the minimum value oft seems 
to be 75. 

We also conjecture that similar properties hold for the generalized q- 
Catalan numbers. 
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